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■ Abstract 
(N ■ 

^ ' In the wake of exploring CP-violation in the decays of B and Be mesons, we perform 

• the straightforward calculation of their nonleptonic decay rates within a relativistic 

, quark model. We confirm that the decays Be DgD^ and Be DgD^ are well 

^SJ \ suited to extract the Cabibbo-Kobayashi-Maskawa angle 7 through the amplitude 

■ relations because their decay widths are the same order of magnitude. In the b — c 
Q . sector the decays B — > DK and Be — > DD lead to squashed triangles which are 

\ therefore not so useful to determine the angle 7 experimentally. We also determine 

Q_i' the rates for other nonleptonic -Bc-decays and compare our results with the results 

^ • of other studies. 
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As was pointed out in [1] and [2,3] the decays D^D^{D^) are well suited 

for an extraction of the CKM angle 7 through amplitude relations. These 
decays are better suited for the extraction of 7 than the similar decays of the 
Bu and Bd mesons because the triangles in latter decays are very squashed. 
The Be meson has been observed by the CDF Collaboration [4] in the decay 
Be J/iplv. One could expect around 5 x lO^*' Be events per year at LHC [5] 
which gives us hope to use the Be decay modes for the studying CP- violation. 

In the case of the Be — > DsD^{D^) decays the relevant amplitude relations 
can be written in the form [2] 
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_ (1) 

where |-D°) = (^|-D°) + |-D°)^ /a/2 is a CP-even eigenstate. The diagrams de- 
scribing the decays -B^ D^D^ and B'^ — > D^D^ are shown in Figs. 1 
and 2, respectively. The color-enhanced amplitude of 5+ — > D^D^ can be 
seen to be proportional to V^^ ■ T^s ~ 0.0029 • exp(i7). At the same time the 
decay amplitude for — > DfD^ is proportional to He • 0.0088 but 

color-suppressed. Simple estimates made in [2] give 

A{B+ ^ D+DO) 

This implies that all sides of the amplitude triangles suggested in [2,6] have 
similar lengths as shown in Fig. 3. It allows one to extract the magnitude of the 
weak CKM-phase 7 from the measurement of the Bf Df + {D^, D^, D^_^) 
decay widths. The method [6] of the extraction of 7 from Eq. (1) is based on 
the parametrization of the amplitudes as 

A{Bt ^ DfW) ^A{B: DJD') = \A\ ^\ 

_ (3) 
^ D-^D^) = 1^1 e'^ e'\ A{B; D^D^)) = \A\ e''^ e'^ , 

where 5 and 5 are the strong final state interaction phases. Introducing the 
notation \A{B+ DtD\)\ = \A+\ and \A{B- D;Dl)\ = \A_\ Eq. (1) 
can be rewritten as 

\A+\^ +\A_\^=\A\^ + \A\^ + 2\A\ \A\ C0S7 cos{6 - S) , 

(4) 

- \A_\^ = 2\A\ \A\ sin J sm{6 - 6) . 
The four solutions for sin 7 are given by [6] 



AiB- 



D-D^ 



A{B- D-DO) 



0{1). 



(2) 



where Y+ 



\A\ + \A\ 



2 \A, 



\A, 



\A\ - \A\ 



Thus, 



P and = 2 

the measurements of the rates of the six decays in Eq. (1) will determine the 
magnitude of 7 with the four-fold ambiguity in Eq. (5). The way to resolve 
the ambiguity was discussed in [6]. 
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In contrast to Be — > DgD, the corresponding ratios for B — > KD and Be 
DD are [2] 



A{B- 



A{B 
MB, 



A{B- 



K-D^ 



A{B+ D+D^) 



A{B- 
A{B; 



A{B- ^ D-D^) 



■■O{0.1), 
C(O.l) 



(6) 
(7) 



which can be seen to lead to squashed triangles which are not very suited to 
measure 7. 

Some estimates of the branching ratios have been obtained before in [7]- [10] 
with widely divergent results. We employ here a relativistic quark model [11] 
to provide an independent evaluation of these branching ratios. 

This model is based on an effective interaction Lagrangian which describes 
the coupling between hadrons and their constituent quarks. For example, the 
coupling of the meson H to its constituent quarks qi and q2 is given by the 
Lagrangian 



J^int{x) ^gHH{x) jdXijdX2FH{x,Xi,X2)q{Xi)V H\Hq{x2) ■ 



(8) 



Here, A^^ and Vh are Gell-Mann and Dirac matrices which entail the flavor 
and spin quantum numbers of the meson field H(x). The shape of the vertex 
function Fh can in principle be found from the Bethe-Salpeter equation as 
was done e.g. in [12]. However, we choose a phcnomcnological approach where 
the vertex functions are modelled by a simple form. The function Fh must be 
invariant under the translation Fh{x + a,xi + a, X2 + a) — Fh{x, xi, X2) and 
should decrease quite rapidly in the Euclidean momentum space. 

In our previous papers [13] we omit a possible dependence of the vertex func- 
tions on external momenta under calculation of the Feynman diagrams. This 
implies a dependence on how loop momenta are routed through the diagram 
at hand. In our last paper [14] and in the present calculation we employ a 
particular form of the vertex function given by 

F ^ \ xf miXi+m2X2\ 

Fh{x, Xi,X2) =6[x ■ ]^h{{xi - X2) ). (9) 

V mi + 1712 / 

where mi and m2 are the constituent quark masses. The vertex function 
Fh{x,Xi,X2) evidently satisfies the above translational invariance condition. 
We are able to make calculations explicitly without any assumptions concern- 
ing the choice of loop momenta. 
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The coupling constants qh in Eq. (8) arc determined by the so cahed compos- 
iteness condition proposed in [15] and extensively used in [16]. The compos- 
iteness condition means that the renormalization constant of the meson field 
is set equal to zero 



n'^(m^) = 



(10) 



where 11^ is the derivative of the meson mass operator. In physical terms 
the compositeness condition means that the meson is composed of a quark 
and antiquark system. For the pseudoscalar and vector mesons treated in this 
paper one has 



np(/) 



■p 



2p2 ^ J A-kH 



tr 

1 
3 

tr 



1 d rd^k ^ 



■P 



2p2 ^ pOi J 47^2^ 



where Wij = mj/{mi + ruj). 

The leptonic decay constant fp is calculated from 



3gp fd^k ~ 



4^ 



1 

my' 



fv ' 



The transition form factors P{pi) P{P2){y{P2)) are calculated from the 
Feynman integrals 



3gpgp> fd k_^^^_^^^^^_^^^y^^^^^_^j^^^^_^^^^2-^ 



An' 



AnH 



•tr 



^9p9v fd^k ~ 



47r2 J 47r2 



:^p{-{k + WisPiY) ^v{-{k + W23P2) 



:i2) 



■tr 



^2(^+ ^2)0'^'5i(^+ ^i)7''53(^)7 • ev 



nip + mv 



A 



where = 7^(1—7°). We use the local quark propagators Si{^) = l/{'mi— }t) 
where rrii is the constituent quark mass. As discussed in [13,14], we assume 
that ruH < mi + 1712 in order to avoid the appearance of imaginary parts in 
the physical amphtudes. This holds true for the hght pseudoscalar mesons but 
is no longer true for the light vector mesons. We shall therefore employ iden- 
tical masses for the pseudoscalar mesons and the vector mesons in our matrix 
element calculations but use physical masses in the phase space calculation. 
This is quite a reliable approximation for the heavy mesons, e.g. D* and B* 
whose masses are almost the same as the D and B, respectively. However, for 
the light mesons this approximation is not so good since the ii'*(892) has a 
mass much larger than the K{494). For this reason we exclude the light vec- 
tor mesons from our considerations. The fit values for the constituent quark 
masses are taken from our papers [13,14] and are given in Eq. (13). 



m„ m, rrif. rriK 

(13) 



0.235 0.333 1.67 5.06 GeV 



We employ a Gaussian for the vertex functuon ^H{k%) = exp(— A;|;/A|^) and 
determine the size parameters A^ by a fit to the experimental data, when 
available, or to lattice simulations for the leptonic decay constants. The nu- 
merical values for Ah are 



A^ Ak Ap Ap, Ab Ag, Age ^^^^ 

1.00 1.60 1.70 1.70 2.00 2.00 2.05 GeV 



We have used the technique outlined in our previous papers [13,14] for the nu- 
merical evaluation of the Feynman integrals in Eqs. (11) and (12). The results 
of our numerical calculations are well represented by the parametrization 

^ ^ 1-as + bs^ ^ ^ 



with s = q^ /m\^. Using such a parametrization facilitates further integrations. 
The values of F(0), a and h are listed in Table 1. The calculated values of the 
leptonic decay constants are given in Eq. (16). They agree with the available 
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experimental data and the results of the lattice simulations. 



fK+ /do /d*o fos 



/d? 



fE 



-(16) 



0.161 



0.215 0.348 0.222 0.329 0.180 0.398 GeV 



The relevant effective Hamiltonian for the decays Be DgD^ and — > DgD^ 
is written as 

Hef£ = {Cl(Ai) iycsV^b • {bu)v-A{cs)v-A + KsKb " {bc)v-A{us)v-A) 

+ C2(/X) {VcsVi ■ ibs)v-Aicu)v-A + K.KI • ibs)v-Aiuc)v-A)} (17) 

where V — A refers to = 7^(1 — 7^). We use the numerical values of the 
Wilson coefficients at the renormalization scale /i = rrib^poie given by Ci = 1.107 
and C2 = —0.248 as given in [17]. Note that we interchange the labeling 1 <-> 2 
of the coefficients to be consistent with the papers [7]-[10]. 

Straightforward calculation of the matrix elements of the decays B^ DgD^iyDgD^ 
by using the effective Hamiltonian in Eq.(17) reproduces the result of the fac- 
torization method. We have 



Gf 

X {a, [f^^^imlj (m|^ - mlo) + /^=^(m?,J ml] ■ fn^ 



(18) 



+ a2 [/f ^^(m^o) {ml - m%) + /5^^(m^o) m],.] • /^o} , 



A(B+ ^ D+DO) = ^ • nlK. (19) 

X a2 [/f ^^(m^o) {m% - m^J + /5^^(m^o) m%] ■ fno 
+ annihilation channel 

where ai = Gi + C,G2 and 02 = C2 + ^Ci with ^ = l/Nc- As usual we put 
the QCD color factor ^ = according to 1/iVc-expansion. Also we drop the 
annihilation processes from the consideration. Note that the calculation of 
the matrix elements of the nonleptonic decays involving the vector D-mesons 
in the final states proceed in a similar way. We extend our analysis to the 
semileptonic and nonleptonic decays of S-meson. 

For numerical evaluation we have used the set of the parameters: mB+ — 5.279 
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GeV, tb+ = 1.655 ps, uib^ — 6.4 GeV, tb^ — 0.46 ps, ai|^=o = 1.107, a2|^=o = 
-0.248 and 



\Vud\ \Vus\ \Vub\ \Vcs\ \Vbc\ 

0.98 0.22 0.003 0.22 0.98 0.040 



(20) 



First, to illustrate the quality of our calculations, we list some branching ratios 
of the i?-meson decays in Table 2 and compare them with the experimental 
data. The exclusive nonleptonic decay widths of the B and Be mesons for 
arbitrary values of ai and 02 are listed in Table 3 whereas their branching 
ratios for ai — 1.107 and 02 = —0.248 are given in Table 4. One can see 
that as it was expected the magnitudes of the branching ratios of the decays 
Be DgD^ and Be — > DgD^ are very close to each other. It gives us hope that 
they can be measured in the forthcoming experiments. Finally, in Table 5 we 
compare our results with the results of other studies. One can see that there 
are quite large differences between the predictions of the different models. 

This work was completed while M.A.I, and O.N. P. visited the University of 
Mainz. M.A.I, appreciates the partial support by GraduiertenkoUeg "Eichthe- 
orien", the Russian Fund of Basic Research grant No. 01-02-17200 and the 
Heisenberg-Landau Program. O.N. P. thanks DAAD and Federal Russian Pro- 
gram "Integracia" Grant No. 3057. 
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Table 1 

Form factors for 5+ — D'^{D*^) and 5+ — > £)+(£)*+) transitions. Form factors are 
approximated by the form F{q^) = F(0)/(1 — as + bs^) with s = /m\^. 





B+ ^ D^{D*^) 


Bt ^ Dt{Dt+) 




m 


a 


b 


m 


a 


b 


F+ 


0.189 


2.47 


1.G2 


0.194 


2.47 


l.Gl 


F_ 


-0.194 


2.43 


1.54 


-0.183 


2.43 


1.53 


Aq 


0.284 


1.30 


0.15 


0.312 


1.40 


0.16 


A+ 


0.158 


2.15 


1.15 


0.168 


2.21 


1.19 


A_ 


-0.328 


2.40 


1.51 


-0.329 


2.41 


1.51 


V 


0.296 


2.40 


1.49 


0.298 


2.41 


1.49 
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Table 2 

Comparison of some branching ratios of the S-meson decays with the available 
experimental data. 





i ri IS \\'< )rK 


\ Dvi! [\-o\ 


7->4. r^fl 4~ 


0.024 


0.0215±0.0022 


B+ D*Oe+u 


0.056 


0.053±0.008 


B+ K+DO 


2.8 • 10"^ 


(2.9 ±0.8) • 10-4 


B+ ^D+W 


0.013 


0.013 ± 0.004 


B+ D+ D*o 


0.008 


0.012 ± 0.005 


B+ D*+D^ 


0.019 


0.009 ± 0.004 


B+ D*+ D*o 


0.046 


0.027 ±0.010 



Table 3 

Exclusive nonleptonic decay widths of the B andSc mesons in 10^^^ GeV. 



B+ ^ K+ L>o 


(0.364 ai ±0.28602)^ 


B+ ^ K+ L>o 


0.00915 


B+ K+ L>*o 


(0.342 04 + 0.442 aa)^ 


B+ K+ 


0.0219 ai 


B+^D+W 


4.367 a? 






B+ D+ D*o 


2.707 a? 






B+ ^ D*+D^ 


6.300 a? 






B+ D*+ D*o 


14.84 a? 






B+ ^ D+ Z)0 


(0.0147ai ±0.014602)^ 


B+ -^D+W 


0.753 ai 


B+ D+ D*o 


(0.0107ai ±0.023402)^ 


B+ D+ L>*o 


1.925 ai 


B+ D*+ 


(0.0233 ai +0.010602)^ 


B+ D*+ W 


0.399 ai 


D+ — D*+ D*^ 


(0.0235 r/i +0.023502)^ 


B+ — D*" L)*" 


1.95 r/i 


Bt^DtD^ 


(0.0689 ai +0.0672 aa)^ 


B+^D+W 


0.0405 ai 


B+ Df 


(0.0503 ai ± 0.106 02)^ 


B+ Df L>*o 


0.101 ai 


B+ D*+ L»o 


(0.101 ai ±0.049802)^ 


B+ D*+ W 


0.0222 ai 


B+ D*+ D*o 


(0.104ai +0.110a2)2 


B+ D*+ L>*o 


0.109 ai 
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Table 4 

Branching ratios of some nonleptonic decay widths of the B &ndBc mesons in cal- 
culated for ai = 1.107 and 02 = -0.248. 



B+ K+D^ 


2.76 • lO"'' 


B+ ^ K+ L»o 


1.41 • 10-6 


B+ - K+ 17*" 


1.82 ■ IQ-^ 


B+ — Iv " 17*" 


3.38 ■ 10-6 


B+ ^D+D^ 


1.11 • 10-'^ 


B+ D+DO 


3.24 • 10-5 


B+ D+ 


0.25 • lO"'^ 


B+ D+ 17*0 


8.28 • 10-5 


B+ D*+ 


3.76 • lO"'^ 


B+ D*+ W 


1.71 • 10-5 


B+ D*+ 


2.84 • lO"'^ 


B+ D*+ 17*0 


8.38 • 10-5 


B+ - Dt 17" 


2. 18 ■ 10"'^ 


B+ - 17+17" 


1.71 ■ 10-6 


B+ D+ D*o 


0.60 • 10-6 


B+ D+ D*o 


4.34 • 10-6 


B+ D*+ 


6.88 • 10-6 


B+ 17*+ 


0.95 • 10-6 


B+ DI+ D*^ 


5.41 • 10-6 


B+ 17*+ 17*0 


4.69 • 10-6 



Table 5 

Exclusive nonleptonic decay widths of the Be meson in units of 10-^5 QeV. 
Comparison with other studies. 



Process 


This paper 


[7] 


[8] 


[9] 


[10] 


[5] 


1?+-17+:D^ 


0.0405 ai 


0.0340 ai 


0.168 ai 


0.01 ai 


0.0415 ai 


0.176 ai 


B+ 17+ D*o 


0.101 ai 


0.0354 ai 


0.143 ai 


0.009 ai 


0.0495 ai 


0.260 ai 


B+ 17*+ 170 


0.0222 ai 


0.0334 ai 


0.0658 ai 


0.087 ai 


0.0201 ai 


0.166 ai 


B+ 17*+ 17*0 


0.109 ai 


0.0564 ai 


0.128 ai 


0.15 ai 


0.0597 ai 


0.951 ai 
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A{B+ ^ DfD^)^AiB- ^ D;D°) 



Fig. 3. The amplitude triangles for the decays {D^ , , D^} . 
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